Anomalous non-ergodic scaling in adiabatic multicritical quantum quenches 
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We investigate non-equilibrium dynamical scaling in adiabatic quench processes across quantum 
multicritical points. Our analysis shows that the resulting power-law scaling depends sensitively on 
the control path, and that anomalous critical exponents may emerge depending on the universality 
class. We argue that the observed anomalous behavior originates in the fact that the dynamical 
excitation process takes place asymmetrically with respect to the static multicritical point, and that 
non-critical energy modes may play a dominant role. As a consequence, dynamical scaling requires 
introducing new non-static exponents. 
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Establishing dynamical scaling relations in many-body 
systems adiabatically driven out-of-equilibrium across a 
quantum phase transition has important implications for 
both condensed-matter physics^, and adiabatic quantum 
computational. A paradigmatic scenario is offered by the 
so-called Kibble-Zurek scaling (KZS}^>^, whereby an ho- 
mogeneous rf-dimensional system is linearly driven with 
a constant speed l/r across an isolated quantum criti- 
cal point (QCP) described by equilibrium critical expo- 
nents v and z. Under the assumption that, in the ther- 
modynamic limit, the system loses adiabaticity through- 
out an "impulse region" [tc — i,tc + i\ which is centered 
around the QCP and has a characteristic width 2t, ex- 
citations are generated in the final state with a density 
nex(t/) ~ T-'^''/(''^+^). While the KZS and its non-linear 
generalizations have been verified in several exactly solv- 
able models^, departures from the KZ prediction may oc- 
cur for more complex quench processes - involving either 
isolated QCPs in disordered^ and infinitely-coordinated 
systems^, or non-isolated QCPs (that is, quantum crit- 
ical regions) ^i^i^° . Evidence of non-KZS, however, has 
also been reported in the apparently simple situation 
of a quench across a single quantum multicritical point 
(MCP) in clean spin chain a^^i^^ . 

In this work, we show how multicritical quantum 
quenches dramatically exemplify the dependence of non- 
equilibrium scaling upon the control path anticipated in 
Ref. i, and demonstrate that anomalous "non-ergodic" 
scaling may emerge in the thermodynamic limit. While 
a non-KZS nex(i/) r~^/^ was previously reportect^i 
and an explanation given in terms of an "effective dy- 
namical critical exponent" Z2 = 3, the meaning of such 
exponent relied on the applicability of a Landau-Zener 
(LZ) treatment, preventing general insight to be gained. 
We argue that the failure of KZS is physically rooted in 
the shift of the center of the impulse region relative to 
the static picture, and that Z2 is determined by the scal- 
ing of a path-dependent minimum gap which need not 
coincide with the critical gap. Furthermore, such a dy- 
namical shift may also cause the contribution from inter- 
mediate non- critical energy states to dominate the scaling 
of the excitation density, via an "effective dimensionality 



exponent" ^2 7^ 0. We show that the latter leads to the 
emergence of new scaling behavior n^^{tf) ~ r~'^/*. A 
unified understanding is obtained by extending the adi- 
abatic renormalization (AR) approach of Ref. |9|. 

Model Hamiltonian.— We focus on the alternating 
spin- 1/2 XY chain described by the Hamiltonian^iii^: 
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where 7± = (l±7)/2, hi ^ h — (— )'<5, and periodic 
boundary conditions are assumed. Here, ft,, i5 e IR are 
the uniform and alternating magnetic field strength, re- 
spectively, whereas 7 G IR is the anisotropy (lifting 
the restriction 7 G [0, 1] is essential for the present 
analysis). An exact solution for the energy spectrum 
of H may be obtained through the steps outlined in 
Ref. 0. The problem maps into a collection of non- 
interacting quasi-particle labelled by momentum modes 
k € K+ = {7r/Af,37r/iV, ...,7r/2 - vr/iV}, whose ex- 
citation gap is given by Afe(7,/i, J) = 4 [ft^ -I- (5^ + 

cos^ fc -I- 7^ sin^ k - 2-y/ h? cos^ k -I- 5'^{h? -I- 7^ sin^ fc)]^/^. 
The quantum phase boundaries are determined by the 
equationa^iii^ h? = 5'^ + \ , 5'^ = It? -\- ^"^ . Thus, the criti- 
cal lines on the 7 = plane consist entirely of MCPs. 

Quench dynamics: Exact results — We assume that 
the system is initially in the ground state, and that (in 
the simplest case) a slow quench across a MCP is im- 
plemented upon changing a single control parameter ac- 
cording to 5\{t) = \{t) - Ac = |(i - ic)/T|"sign(t - tc) 
over a time interval t £ [tQ,tf], where a — 1 corresponds 
to a linear driving, and Ac is the critical value. Thus, 
the time-dependent Hamiltonian H{t) may be written as: 
H{t) ~ Hc + SX{t)Hi, where He is quantum-multicritical 
at time t^ in the thermodynamic limit, and Hi is the con- 
tribution that couples to the external control (a similar 
parametrization is possible for quenches involving mul- 
tiple parameters). Without loss of generality, we may 
let tc = 0. In what follows, we shall focus on two rep- 
resentative MCPs, A and B as marked in Fig. [TJ each 
approached through two different paths (path 5 will be 
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introduced later) with the following properties: 
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FIG. 1: Phase diagram of H in Eq. ([T]) when h — 1 (left) and 
5 = (right). The dashed-dotted (green) line separates the 
ferromagnetic (FM) and paramagnetic (PM) phases, the solid 
(red) lines separate dimer (DM) and FM, whereas the dotted 
(green) line is the superfluid phase (SF). The arrows indicate 
the control paths we choose to approach the MCPs A and B. 

In order to quantify the amount of excitation at a 
generic instant t, we numerically integrate the time- 
dependent Schrodinger equation for H(t) and monitor 
two standard "non-adiabaticity" indicators^i^ii^: the ex- 
citation density, n^^, and the residual energy, AH. For 
a linear quench along either path 1 or 2 (left panel of 
Fig. El), we find that n^^{t) - r-'^/(''^+i) = t^^/^ and 
AH{t) ~ = r-\ which is consistent 

with KZS^ and our conclusion in Ref. HI. For paths 
3 and 4, however (right panel of Fig. [2]), we find that 
nex(ti~ T-i/*^ and AH{t) - t'^/s, which is non-KZS (in 
Ref. M, the t^^/^ scaling was pointed out for an equiv- 
alent quench scheme across MCP A). Similar anomalous 
exponents are found for non-linear quenches along paths 
3 or 4, e.g., ne^{t) ~ r~^/^ for a = 2. 

The above results clearly show that, for quenches 
across a MCP, whether KZS is obeyed depends sensitively 
on which control path is chosen. A closer inspection re- 
veals the following important differences: (i) Paths 1, 2 
start and end in essentially the same phase, correspond- 
ingly the excitation spectrum is invariant under a trans- 
formation A 1-^ —A of the control parameters. Paths 3,4 
do not exhibit this symmetry; (ii) Along paths 3, 4, the 
MCPs A and B belong to the Lifshitz universality class 
(i^ = 1/2), although all paths share z = 2. It is then 
natural to ask which of these differences may play a role 
in determining the anomalous dynamical scaling behav- 
ior. To answer this question, we introduce another path 
across MCP A (path 5), h{t) = l + |7(t)| = l + \t/T\, which 
starts and ends in the PM phase but, in each of the two 
segments, crosses the MCP A with Lifshitz exponents. 
Surprisingly, the observed scaling is nex(i) ~ t~^^^ (left 
panel of Fig. [3]), which is neither KZS nor —1/6. An 
identical —3/4 scaling holds for a "V-shaped" path across 




FIG. 2: Exact scaling of the excitation density throughout a 
linear quench along path 2 (left) and path 3 (right). Right 
inset: Scaling of the residual energy along path 3. 




FIG. 3: Left: Exact scaling of the excitation density through- 
out a linear quench along path 5. Right: Scaling of the final 
excitation density in path 5 for different system size: nex(t/) is 
the same to a numerical accuracy of 10"*^ , up to r < 2 x 10^ . A 
linear fit yields -0.747± 0.001 over the range 200 < r < 2000. 

MCP B, that starts and ends in the DM phase. As finite- 
size analysis reveals, all the observed anomalous scalings 
are practically independent upon system size over a wide 
range of quench rates (see e.g. right panel of Fig. [3]), 
establishing them as truly thermodynamic in natureii. 

Landau-Zener analysis.— We begin to seek an under- 
standing from limiting cases where an exact solution 
for ne^{tf) may be obtained based on the LZ picture^^. 
This is possible provided that a = 1 and the Hamilto- 
nian can be decoupled into effective two- level systems. 
Among the above-mentioned paths, only paths 4 and 5 
(for which (5 = 0) can be exactly mapped to a LZ prob- 
lem, thanks to the possibility of rewriting H in Eq. |(T|) as 
H = EkHk^ Efc BlHkBk, where bI = (c_fe, 4) and 

TT _(Hk.ii Hk,i2\_ f-h + cosk 7sinfc \ , , 
\H*k.i2 -HkMj^ 7sinfc h-cosk)-^^> 

A rotation Rk{qk), Qk € [— 7r/2, 7r/2), renders the off- 
diagonal terms in Eq. ([2]) independent upon 7 (hence 
time), allowing use of the LZ formula. Consider path 
4 first. By choosing tan2gfe = — sinfc, the transformed 
Hamiltonian matrix elements become H'^. ~ — 2(1 — 
cosfc) cos2qfe — 2i/r(cos2g/j — sinfcsin2gfc), and H'j^ = 
2(1 — cosfc) sin 2gfe. If the critical mode fcc is defined by 
requiring A^.^ = in the thermodynamic limit, we have 
fcc = for the MCP A. We may then let tan 2qk ~ sin 2^^, 
and the appropriate qk ~ — fc/2. From the LZ formula, 
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the asymptotic {tf oo) excitation probability reads 

— 27rr (1 — cos k)^ sin^ 2qt / (cos 2qt —sin k sin 2gfc ) ^ —7TTk^/2 

Pk ^ ~ ^ 7 

where the approximation follows from a Taylor expansion 
around kc. Integrating over all modes yields nf.y^{tf) 
T~^l^ , which is consistent with our exact numerical re- 
sult. Therefore, mathematically, the t~^I^ scaling follows 
from the fact that the exponent in scales as k^ — k'^^^ , 
with Z2 = 3. In turn, this originates from the scaling 
of the off-diagonal terms Hj^ -^^ ^ • Physically, as we 
shall later see by invoking AR, H'^ m^iy be interpreted 
as the minimum gap for mode k along path 4. 

To unveil the t~^/^ scaling, it is necessary to use the 
exact finite-time LZ solution. For simplicity, we restrict 
to half of path 5, by quenching the system from the PM 
phase up to the MCP A. This has the benefit of avoiding 
the non-analytic time-dependence of the control param- 
eters that path 5 exhibits at A, while leaving the t~^/^ 
scaling unchanged thanks to the symmetry of the excita- 
tion spectrum. Starting from Vitanov's expression [Eq. 
(7) in Ref. the excitation probability Pk{tf ) can be 
computed via the parabolic cylinder function D^{z), 



Pkitf) 



D,^y2{TfV2e''''^') cos 0{Tf) 
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where w = (1— cos k) sin 2qky/T/ -y/cos 2qi; + sin 2qk sin k ^ 
k^y/r is the rescaled coupling strength, Tf = —Loj sinfc ~ 
— k'^y/r is the rescaled time, tan2gfc — sinfc, and 
0{Tf) = 1/2 arctan (tj/T/) -I- 7r/2. Since for our quench 
process a; <C \Tf\ ^ 1 around fcc, we may estimate Pk{tf) 
by Taylor-expanding Dy{z) around T/ = 0: 

Pk{tf) « (l-e-""'/2)/2 + cos2 0(T/)e-""'/2 

- sin20(T/)/2sinxfcVl-e---", (3) 

where Xfc ~ 7''/4 around kc- This approximation 
breaks when Tf ^ 1, setting the scaling of the highest- 
momentum contributing mode, that is, fcmax ~ r~^/^. 
In Eq. the dominant term cos^ 6'(T/)e~'''^'/2 



fc2 since e'""^'/^ 



1 within 



^max, which means Pk{tf) k^. Thus, nej^{tf) = 

jQ""^^Pk{tf) ^ ^max ~ T~^^^, in agreement with our 
numerical results. Remarkably, the fact that Pk{tj) ^ 
{k—kcY^, d2 = 2, indicates that fcc is noi excited despite 
a static QCP being crossed, and also that the excitation 
is dominated by intermediate energy states. In fact, at 
the MCP A, the modes around fcc are still far from the 
impulse region, since |T/| uj, which sets the LZ tran- 
sition time scalei^. This is in stark contrast with the 
main assumption underlying KZS, where the center of 
the impulse region is the static QCP, and excitations are 
dominated by modes near fcc, as reflected in the typical 
scaling pfc ~ (fc — fcc)°. Therefore, the shift of the actual 
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FIG. 4: Scaling of the excitation density from first-order AR 
for a linear quench along path 3 (left) and half-5 starting at 
MCP A (right). Right inset: Low-lying single-mode excitation 
spectrum along path 4 for A'' = 100. 



(dynamical) impulse region relative to the static one is 
ultimately responsible for the anomalous r~'^/^ scaling. 

Perturhative scaling approach.— Since the system be- 
comes gapless at a single MCP along all the paths un- 
der study, first-order AR is a viable approach^. Let 
iV'm(O) be a basis of snapshot eigenstates of Hit), with 
snapshot eigenvalues Em{t), m = labelling the ground 
state. The time-evolved state may be expanded as 

m)) = (^o\t)\Mt)) + E™^oCm Wl^mW), where the 
coefficients Cm\t) determine the excitation amplitudes 
and are given by Eq. (4) in Ref. i. First-order AR cal- 
culations of nex(i) demonstrate that for linear quenches 
along paths 1 and 2, riex ^ t~^^^, whereas riex ^ t~^I^ 
along paths 3 and 4 (left panel of Fig.|4]). Since the non- 
analyticity at A in path 5 might cause problems in AR, 
again we choose to study half of path 5 (right panel of 
Fig. H]) . All the AR results agree with the exact simula- 
tion results, confirming that AR reproduces the correct 
dynamical scaling across a generic isolated QCP. 

Predicting the scaling exponent based on AR requires 
scaling assumptions for the contributions entering Cm^ (t) 
{i.e., A™(t) = Er,^{t) - Eo{t) and {i;„,{t)\Hi\Mt))] and 
the ability to change discrete sums of all the contribut- 
ing excited states into integrals, for which the density of 
excited states p{E) is required. Since typically the AR 
prediction is consistent with KZS, anomalous behavior 
must stem from anomalous scaling assumptions of (one or 
more of) these ingredients. We first examine the excita- 
tion spectrum along different paths. Notice that since Hi 
is a one-body perturbation, only single-mode excitations 
are relevant, thus the index m labeling many-body exci- 
tations may be identified with a momentum mode. Along 
paths 1 and 2, it turns out that the minimum gap among 
all modes is always located at fcc, whereas along paths 3 
and 4, the minimum gap is located at fcc only at the MCP. 
This suggests that knowing the critical exponents of the 
MCP alone need not suffice to determine the dynami- 
cal scaling due to the existence of "quasi-critical" modes 
along paths 3 and 4. In fact, mathematically, along path 
4, 9Afc(7, 1 -I- 7, 0)/dj = gives the location of the min- 
imum gap for each mode fc at 7 = (cos fc— 1) /(I -I- sin fc^), 
which is largely shifted into the FM phase (see inset in 
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Fig. 131). By inserting this relation back into Afc, the 
composite function Afe(7) = Ak ^ {k — kcY . Fol- 
lowing the same procedure also yields Kk ^ {k — kcY 
along path 3, while has the same scaling as A^ at 
the MCP along paths 1 and 2. This motivates modi- 
fying the AR scaling assumptions of Ref. as follows: 
E^it) - Eo{t) = SXitr^ f„,{A,n{tmin)/SX{tr'), whcrc 
Am(imin) IS the minimum gap of mode m attained at 
tmin along the path, and fm is a scaling function. 

The above modification requires the scaling of p{E) to 
be modified by letting p{E) ^ E'^^^^~^^ where Z2 comes 
from the dispersion relation of Am(tnun)- If the mini- 
mum gap of any mode is below a certain energy along the 
path, that mode should be counted into the contributing 
excited states. Accordingly, we have Z2 — z ~ 2 along 
paths 1, 2, half-5, and Z2 = 3 along paths 3, 4. Back to 
the LZ analysis, note that the off-diagonal term H[2{k) 
is the minimum gap of mode k along the path if there 
exists a time at which the diagonal term H{i{k) = 0, 
as it happens for path 4. For path 5, however, the off- 
diagonal term never becomes the minimum gap since the 
system never leaves the PM phase. Therefore, the off- 
diagonal term in the LZ picture need not suffice to de- 
termine the dynamical scaling, and the shift of the loca- 
tion of the minimum gap for each mode from the static 
QCP is at the root of the anomalous behavior we ob- 
serve. Lastly, we consider the matrix elements of Hi. 
Numerical simulations suggest that {il'm(t)\Hi\ipQ{t)) — 
6X{ty^-^gm{Am{t^in)/6X{ty), where Qm is a scaling 
function, and Am(tmin) is the minimum gap of mode m 
along a path that extends the actual path to i/ oo 
when the quench is stopped at the MCP, and coincides 
with the actual path otherwise. Then along paths 1 and 
2, Am(tniin) ^ k^ , whcrcas along paths 3, 4, and half-5, 
Am(imin) ~ k^ . Together with the other scaling assump- 



tions, and taking the linear case a = 1 as an example, AR 
yields |c^^| - k°, n^^ - ^-iz/z2){,^/i,^z+i)) ^^^^^^ p^^^j^g ^ 
to 4, and Ic^^l - fc\ n^^ - ^-3^/(<..+i) g^i^^^ ^g^^^i, 
which completely agrees with the numerical results. 

Building on the above analysis, we argue on physi- 
cal grounds that the scaling of the excitation density for 
quenches across an arbitrary (standard or multicritical) 
isolated QCP is determined by three conditions: (i) From 
the condition of adiabaticity breaking, the typical gap, 
A, scales as A ~ ^-auz/{ai^z+i).^ ^^-j accessible ex- 
cited state contributes to the excitation if and only if its 
minirnum gap along the path matches with this typical 
gap, Ak ~ A, with Afc ^ (/c — kcY^\ (iii) The excitation 
probability, pk, scales as Pk ^ {k — kcY^'^ , where c?2 can 
differ from if the center of the impulse region is greatly 
shifted relative to the static limit. Then upon integrating 
up to energy A, and using pe ^ Pk{E) ^ E'^^^^^ , yields 

riex ^{d,+d2)/z2 ^ ^-{d+d2)ai'z/[z'2(avz+l)] ^ 

which is consistent with all the results found thus far. 
While KZS corresponds to Z2 = z, c?2 = 0, situations 
where Z2 ^ z and/or d2 ^ Q are genuinely dynamical: 
knowledge about the path-dependent excitation process 
becomes crucial and non-equilibrium exponents cannot 
be fully predicted from equilibrium ones. Interestingly, 
in the model under examination the Lifshitz universality 
class appears to be the only universality class for which 
anomalous scaling occurs, among all possible paths in- 
volving MCPs. Whether Lifshitz behavior may consti- 
tute a sufficient condition for anomalous behavior re- 
quires further investigation in other many-body systems. 
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